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Measuring the equation of state of trapped ultracold bosonic systems in an optical 

lattice with in-situ density imaging 

Ping Nang Ma/ Lode Pollet,^ and Matthias Troyer^ 

^ Theoretische Physik, ETH Zurich, 8093 Zurich, Switzerland 
^Physics Department, Harvard University, Cambridge 02138, Massachusetts, USA 

We analyze quantitatively how imaging techniques with single-site resolution allow to measure 
thermodynamical properties that cannot be inferred from time-of-light images for the trapped Bose- 
Hubbard model. If the normal state extends over a sufficiently large range, the chemical potential 
and the temperature can be extracted from a single shot, provided the sample is in thermodynamic 
equilibrium. When the normal state is too narrow, temperature is low but can still be extracted 
using the fluctuation-dissipation theorem over the entire trap range as long as the local density ap- 
proximation remains valid, as was recently suggested by Qi Zhou and Tin-Lun Ho [arXiv:0908.3015]. 
However, for typical present-day experiments, the number of samples needed is of the order of 1000 
in order to get the temperature at least 10% accurate, but it is possible to reduce the variance by 
2 orders of magnitude if the density-density correlation length is short, which is the case for the 
Bose-Hubbard model. Our results provide further evidence that cold gases in an optical lattices can 
be viewed as quantum analog computers. 



I. INTRODUCTION 

The prototypical, unsolvable models of condensed mat- 
ter physics can be implemented in a tunable, clean and 
controllable way with cold atoms in optical lattices [l[. 
Experiments have the potential to offer new insight in 
the long-standing problems of condensed matter physics, 
and would ultimately allow us to establish the validity of 
those models in the description of more complicated ma- 
terials. Before such quantum simulators can be trusted 
they need to be validated by benchmarking them against 
known results of models that can be solved accurately on 
a classical computer. This was done in great detail for 
the superfluid to normal liquid transition at unity filling 
in the Bose-Hubbard model by comparing experimentally 
observed interference patterns to the ones computed in 
full ab-initio quantum Monte Carlo simulations |la |. 

One important experimental issue is the accurate de- 
termination of temperature. A thermometer measures a 
property of a subsystem in thermal equilibrium with the 
rest of the system. The measurement of this property can 
be used to determine the temperature if the temperature 
dependence of this property is known from theoretical 
calculations or experimental calibration. The accuracy 
of the thermometer increases and the number of mea- 
surements necessary decreases if more information about 
the system is available. 

In this context, recent experimental progress in single- 
site resolution and addressability [9l4l3j provides us with 
important additional tools to validate the Bose-Hubbard 
model. With this it is possible to determine temperature 
in the lattice system using the fluctuation-dissipation 
theorem without invoking theory |19l - [2l| . Also the equa- 
tion of state and the chemical potential can be extracted 
from the edges using high temperature series expansions, 
provided those regions are large enough. This was pre- 
viously suggested when the edges are in the ideal gas 
regime [3, Q; however this is too restrictive. We note 



that number fluctuations were previously suggested as 
an effective thermometer [IJ, [l^. A failure to extract 
temperature and/or chemical potential signals that the 
experiment is not in equilibrium. 

In this paper, we extend our ab-initio study of ultra- 
cold gases [131 to the physics accessible through single- 
site resolution detection tools. We focus on extracting 
the temperature from the lattice experiments. In doing 
so, we extend the studies initiated in Refs. [2|, Q from the 
ideal gas to the strongly interacting systems for realistic 
system sizes, and hereby validate and improve on their 
proposals where applicable. Carefully choosing the range 
over which correlations are measured wc show that accu- 
rate thermometry can be performed with a few dozen 
measurements. Using additional theoretical input, such 
as a second order high temperature expansion of the den- 
sity, it becomes possible - in some parameter regimes - to 
extract the temperature accuractely from a single density 
profile. 

The structure of this paper is as follows. First, we in- 
troduce the model in the next section. Second, we discuss 
a general thermometry scheme based on the fluctuation- 
dissipation (FD) theorem, and illustrate its validity for 
strongly correlated bosons in an optical lattice. Third, 
we investigate how reliably 2™^ order high temperature 
series expansions (HTE2) can be applied to the edges in 
order to extract the temperature and the chemical po- 
tential, before concluding in Section IV. 



II. ULTRACOLD BOSONIC OPTICAL LATTICE 

Strongly interacting bosons harmonically trapped in 
an isotropic optical lattice can be quantitatively de- 
scribed by the single-band boson Hubbard Hamilto- 
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(1) 
where b (6^) and h are annihilation (creation) and number 

operators, respectively. The operator N counts the to- 
tal number of particles found for a system with chemical 
potential /x, and the external harmonic trapping is char- 
acterized by Vt- The nearest neighbor hopping t, the on- 
site repulsion strength U, and the confinement strength 
Vt are derived from the lattice laser potential Vq through 
band structure calculations j22| . 

To mimic in-situ density measurements in ultracold 
bosonic optical lattice experiments, we have performed 
Quantum Monte Carlo worm simulations [23, [23 (exact 
up to the statistical errorbars) and obtained a time- 
series of 3-dimensional density measurements which are 
then column-integrated along the line-of-sight direction. 
Uncorrelated measurements arc obtained by imposing a 
strict criterion for the autocorrelation time r < 0.2 on our 
simulation data used in subsequent thermometry analy- 
sis. 

Throughout this paper, we focus on a physical system 
of 125,000 ^^Rb atoms. We choose realistic parameters, 
using a trapping frequency Vr/t = 0.0091 for U/t = 10 
and Vr/t = 0.0277 for U/t ~ 50. For convenience, we 
use the lattice spacing a = A/2 as unit of length. 



III. THERMOMETRY SCHEME BASED ON 
THE FLUCTUATION-DISSIPATION THEOREM 

In this section, we illustrate how the fluctuation- 
dissipation (FD) theorem can be turned into an effective 
temperature probe for ultracold bosonic optical lattices, 
enabled by in-situ density imaging experiments [9|, Il0| . 
Our ap pro ach is a generalization of the proposal of Zhou 
and Hofl^. It is based upon the 3D density. 



(n(r)) = {n{p,(j),z;T,^j.)) = 



Trn(p,(/),z)e-'^(-^-^^) 
Tre-'3(^-A'^) 



(2) 

integrated along the line-of-sight {n{p)) = Jdz{n{r)), 
and the integrated density-density correlations, 

RdP) -j'^P'i {n[p)n{p'))-{n{p)){n{p')) } 9 {£, - \p - p'\ 

(3) 
within a window size ^. Both quantities can be mea- 
sured directly from in-situ density images. Here, r de- 
notes a coordinate in three dimensions parametrized in 
cylindrical coordinates as r(p, 0, z), while p(p, 0) denotes 
the in-plane coordinate. 6'(.) is the Heaviside step func- 
tion. Under the assumptions of the validity of the lo- 
cal density approximation (LDA) and a value of ^ that 
is larger than the density-density correlation length, the 
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FIG. 1. Illustration of the quantities entering the fluctuation- 
dissipation thermometry formula. Shown from top to bot- 
tom are ; 1) cross-section density nir) , 2) column integrated 
density n{p) , 3) dissipation term L{p) , 4) fluctuation term 
(5 — 3) Raip) , 5) fluctuation term (^ — oo) Roo{p)- We take a 
3D bosonic *^Rb optical lattice system with A'^ = 125, 000 and 
we average over 1000 independent measurements obtained 
from a QMC simulation. The parameters in the left column 
are U/t = 10 , T/t = 1; and in the right column: U/t = 50 , 
T/t = 1. 



FD theorem takes the form T x L{p) = R^{p) where the 
I (iissipation term 



L{p) 
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dp 
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can be computed in the way shown in Appendix]^ The 
LDA is a very good approximation for the density pro- 
file and only breaks down in the vicinity of the critical 
point [13, [i3- After averaging over the angular vari- 
able (j) for radially symmetric lattices, ie. {Ljp),Rjp)} = 
^ J d4>{Ljp),Rjp)}, the temperature T can be esti- 
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mated from a least-square fit taking the measurement 
errors into account: 
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FIG. 2. Illustration of the fluctuation-dissipation thermom- 
etry scheme by showing difl^erent L{p)-R^{p) plots at various 
window sizes ^ = 0, 1, 2, 3, 5, cxd. This approach illustrates 
how the density-density correlation length can be found in 
an experimental system. When the window size is smaller 
than density-density correlation length, systematic errors set 
in and this results in nonlinear L{p)-R{p) behavior, while for 
^ larger than the density-density correlation length the be- 
havior of L{p) — R{p) is linear. However, statistical noise also 
increases with increasing window size. We take a 3D bosonic 
*''Rb optical lattice system with N = 125, 000 and we average 
over 1000 independent measurements obtained from a QMC 
simulation. The parameters in the left column are U/t = 10, 
T/i = 1; and in the right column U/t = 50, T/t = 1. 




(5) 



where A^ and An are the errors in L(p) and R{p), re- 
spectively. Here, p is discretized into bins with width 
corresponding to the experimental resolution. State of 
the art experiments allow us to set the binwidth to unity, 
even though we find that the scheme can tolerate a res- 
olution of up to 5 sites (see Appendix |B]) |9l-[l3l| . In the 
limit S. — !■ oo, this scheme reduces to the one of Zhou and 
Ho [l9| in which the authors obtained an estimate for the 
temperature in an optical lattice of 1200 no n- interacting 
fermions with 3% error over 50 independent samples for 
a temperature T/t = 0.1. 

For fermions, both Ref. [201 ^''^^ ^^^- 1211 suggested 
that the FD theorem is an absolute thermometer in an 
harmonic trap without lattice, through careful normal- 
ization [201 of Rooip), However, for the small number 
0(10^) — O(IO^) of fermions in their experiments, the 
estimated temperature deviated from the time-of-flight 
(TOF) measurements by about 30% [2l|. 

For bosonic optical lattices, our scheme is illustrated 
in figure [Tj and [2] for the temperature of T/i = 1, whereby 
a window size of ^ = 3 is sufhcent to capture (almost) 
all the correlations. Using ^ = 3 as the standard for 
T/t = 1 and higher temperatures, we show the number of 
independent measurements needed to estimate the tem- 
perature within 5% error in Table HI The enormous vari- 



System 


nr of shots 


C=3 


C=oo 


U/t = 10 , T/i = 1 
U/t = 10 , T/i = 3 
U/t = 50 , T/i = 1 
U/t = 50 , T/t = 3 


20 
14 
21 
12 


0(10*) 
0(10*) 
0(10*) 
0(10*) 



TABLE I. Number of uncorrelated shots needed to obtain 5% 
error in 3D *^Rb optical lattice experiments trapping 125,000 
bosons with bin-width = 1.0 at U/t = 10 and U/t = 50. The 
variance reduction through window-sizing leads to orders of 
magnitudes improvement. 

ance reduction through window-sizing turns the FD ther- 
mometry scheme into a feasible tool for ultracold bosonic 
optical lattice experiments. With 20 independent shots 
uniformly distributed over 20% spread in T and 1% in 
A^, this scheme remains applicable (see Appendix [C]) . At 
lower temperature, the scheme remains valid although a 
larger ^ is needed, and the statistical noise will inevitably 
grow. An example of T/t = 0.5 is shown in figure [3] where 
the correlations are effectively captured by a window size 
of ^ = 5, and 100 independent measurements are required 
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FIG. 3. Fluctuation-dissipation thermometry scheme at 
shghtly lower temperature. We take a 3D bosonic *^Rb opti- 
cal lattice system with parameters A'' — 125, 000, U/t = 10, 
T/t = 0.5, (, = 5, and we average over 100 independent mea- 
surements obtained from a QMC simulation. 



to attain an accurate temperature estimate with 5% er- 
ror. 

Current optical lattice experiments using fluorescence 
techniques can only measure the parity (even/odd) of the 
occupation number per site. This affects the FD ther- 
mometry scheme, as illustrated in figure |4l but through 
selection of those points in the L{p)-R^{p) that are on 
the linear slope, we could still obtain an acceptable esti- 
mate of the temperature. Deep in the edges, the number 
of doublons is very low compared to the number of holes 
due to the low overall density and the high potential en- 
ergy cost of creating a doublon. 



System 


estimated temperature [t] 


U/t = 10 , T/t = 1 
U/t = 50 , T/t = 1 


0.985 ± 0.008 
1.003 ± 0.012 



TABLE II. The estimated temperature for a 3D bosonic Rb 
optical lattice system in the presence of doublon-hole indis- 
tinguishability. The parameters are A^ = 125, 000, bin-with = 
1.0, 5 = 3 at C//t = 10,50 and T/t = 1. We average over 1000 
uncorrelated measurements obtained from a QMC simulation. 



IV. SECOND-ORDER HIGH TEMPERATURE 
EXPANSION 



U/t= 10,77f=1.0 

.a.8^„ - distinguishable doublons 

°°'°,^ " indistinguishable doublons - 



10 20 30 40 50 60 



m 

60 


'"•■•-.. 










- 




40 


'="""..»„,^ 










- 




20 






°°a 






- 




10 


20 


30 
P 


40 


50 


60 




i(p) 



1.5- 

S i"« 

0.5- 



^ = 50, 77r=1.0 



= dislingushable doublons 
= indistinguishable doublons 



o,t- 



10 20 30 40 50 60 

70| 

60 "■■■■....._ 
50 - 



~40- 
o. 

*30- 
20- 
10- 



10 20 30 40 50 60 

P 

2| 



1.5- 
0.5 




i.(P) 



FIG. 4. Fluctuation-dissipation thermometry scheme in the 
presence of doublon-hole indistinguishability. From top to 
bottom are shown : 1) cross-section density n(r) , 2) col- 
umn integrated density n(p) , 3) dissipation term L(p) , 4) 
fluctuation term (^ = 3) Raip) , 5) L{p)-R3,{p) relation- 
ship. We take a 3D bosonic *^Rb optical lattice system with 
A^ = 125, 000, and we average over 1000 independent measure- 
ments obtained from a QMC simulation. In the left column 
the parameters are U/t — 10 , T/t = 1; and in the right col- 
umn U/t — 50 , T/t = 1. Blue circles (green squares) show 
the curve where doublons can (cannot) be distinguished from 
holes. 



In this section we propose an alternative scheme for 
thermometry based on second-order high-temperature se- 
ries expansion (HTE2). Deep enough in the edges, there 
will always be a normal region for any temperature and 
interaction strength where the system is well described 
by second order high-temperature series expansions. One 
advantage of this scheme is that it works not only for 
integrated column densities but also for density profiles 
measured only in 2D cross sections of a 3D system [? 



]. In addition this scheme allows determination of the 
chemical potential. 

In practice, a single shot will suffice to extract the 
quantities of interest. When HTE2 applies, it also gives 
a foundation for the limit on the small window size £, for 
the former fluctuation dissipation thermometry scheme 
(indeed, there is no point in applying the FD scheme 
whenever HTE2 works) . It is possible to use higher-order 
schemes than HTE2, but we find that the gain is minimal 



compared to the additional effort. However, it may be 
that the density in the edges is so low that it cannot be 
measured because of the low signal-to-noise ratio, which 
will happen for very low temperatures. In this section, 
we wish to make these ideas more quantitative. 

In the absence of correlations the density {ni ) is cap- 
tured by the zcroth-order high temperature expansion 
theory and given by 



(-r^) 



z^ 






where the zeroth order partition function Z^ is 






(6) 



(7) 



and the onsite diagonal energy Di is 

D, = —ni{ni - 1) . 
The zeroth order expansion, suggested in Ref. [21, does 
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~ order (true) 
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FIG. 5. Illustrating the concept of wing thermometry for 
a 3D bosonic *^Rb optical lattice system, le., describing the 
normal state by high temperature series expansions. Blue 
circles: In-situ density profile obtained from 100 uncorrelated 
measurements obtained by a QMC simulation with param- 
eters U/t = 10 , T/t ^ 3 , N = 125, 000. The superfluid- 
normal phase boundary occurs at the density (n) = 0.42 or 
chemical potential fi/t = —2.75. The 2"'^ order series cap- 
tures all the physics in the normal regime, whereas the O"* 
order has a very small validity range. 



(+■-) 

E 



nf7i,e-^(^'+'°^-^-"--''^' 






(10) 



where {i,j) denotes the sum over nearest neighbors. 



(— ,+) stands for (n* ~ Ui — l,n* 
versa for (+,—). F*- is defined as 
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and we also introduce x^j^ 
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with7]J'' = U{n\"'—n^j ' + l) + {fii—fij). The temperature 
is found by fitting the cross-sectional experiment in-situ 
density measurement against the density found in second 
order high temperature series expansions (HTE2) , 



{n^) = (nf ^) + 



^ 7(0)7(0) 



(13) 
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(0)\ 



HTE2 gives a fairly accurate 



description of the edges in present experiments. In prac- 
tice, one should fit from some distance ri till the end 
of the trap, and vary ri in order to find the range of 
applicability of the HTE2 scheme, see Fig. [6] 

If the entire regime of the bosonic optical lattice is in 
the normal phase, one would require no more than a sin- 
gle experimental density measurement to extract a reli- 
able estimate of the temperature and chemical potential. 
This is shown in figure [71 



a rather poor job in describing the edges of the system, 
as illustrated in Fig. \5\ However, sufficient accuracy over 
a wide density range in present-day experiments is found 
by emplying the second order (up to (/3t)^) partition 
function. 
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V. CONCLUSIONS AND 
ACKNOWLEDGEMENTS 

We have discussed how single-site resolution detection 
tools can be used to obtain the equation of state and/or 
the temperature in trapped ultracold gases in an optical 
lattice. We analyzed whether the fluctuation-dissipation 
theorem can be used to extract temperature for present 
experiments when the LDA is valid [T^ over the en- 
tire tap. Taking advantage of the fact that the density- 
density correlation length is short away from the critical 
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FIG. 6. Second order high temperature series expansion ther- 
mometry scheme for a bosonic optical lattice system. 100 
cross-sectional density measurements are used to estimate 
the temperature and chemical potential. The system con- 
sists of a 3D optical lattice with *^Rb atoms with parameters 
iV = 125, 000, At/i = 4.835, U/t = 10, T/t = 3. The blue cir- 
cles are data averaged over 100 measurements obtained from 
a QMC simulation; the red line is a least-square fit over the 
normal region where nat/t = 5.244 and Tat/t = 2.820. 
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FIG. 7. Second order high temperature series expansion 
thermometry scheme for a bosonic system that is entirely 
in the normal phase. No nrore than a single shot of cross- 
sectional density is needed to estimate the temperature and 
chemical potential within 10% accuracy. We take a 3D opti- 
cal lattice system with bosonic ^'^Rb atoms and paranreters 
iV = 125, 000, fi/t = 25.97, U/t = 50, T/t = 3. The blue 
circles are obtained from a single measurement in a QMC 
simulation; and the red line is a least-square fit over the en- 
tire normal region where nat/t = 25.92 and T^t/t — 2.824 
nK. 



region (also in the superfluid phase), a few dozen mea- 
surements arc sufficient in order to extract the tempera- 
ture reliably and accurately under the condition that the 
shape of the chemical potential landscape is known, the 
system is in thermodynamic equilibrium, and that the 
local density approximation holds. 

Using more theoretical input, such as density profiles 
obtained in a second order high temperature expansion 
the temperature can be obtained from the normal edges 
sometimes already with a single measurement. In cases 
where the nornral region on the edge is too narrow one 



can either go to higher order in the high temperature 
expansion or experimentally shape the trap to obtain a 
wider normal region. 
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for parallelization and Monte Carlo data analysis. The 
simulations were performed on the Brutus cluster of ETH 
Ziirich, and the entire workflow was carried out in the 
Vistrails framework. We thank M. Cheneau, T-L. Ho, N. 
V. Prokof'ev, B. V. Svistunov, L. Tarruel, D-W. Wang, 
aird Q. Zhou for useful discussions. We acknowledge fi- 
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Appendix A: Numerical approximation to the 
density derivative 

In the dissipation term 



1 1 djnjp)) 
HP) = -7^77 7{ 

2Vt P op 



(Al) 



the derivative needs to be taken numerically, which al- 
ways involves an approximation. The most direct way 
is to use central differences, but to improve on the qual- 
ity of data one could do the following: first perform a 
local parabolic fitting on each {n{p)) density bin over a 
fitting range of (2k-f 1) density bins, ze., over the interval 
[{n{p — fc)}, {n{p -f k))] before taking the spatial deriva- 
tive analytically from the fitted density profile. In prac- 
tice, we find that the fiuctuation dissipation thermomtery 
scheme is unaffected by choosing different numerical ap- 
proximations to calculating the derivative. However, 
large density gradients always result in systematic errors. 



Appendix B: In-situ density imaging with few-sites 
resolution 



Although we aim at single-site resolution detection 
tools in the analysis of the fluctuation-dissipation ther- 
mometry scheme in this paper, the scheme remains appli- 
cable when the resolution is just a few sites, as illustrated 
in Fig. [5] Here, we simulate such experiments with dif- 
ferent resolutions up to 5 sites, by varying the bin. The 
fiuctuation dissipation thermometry scheme suffers from 
increasing systematic errors in the dissipation ternr L{p) 
both for U/t = 10 and U/t = 50 at T/t ^ 1. Yet, 
a relatively linear L{p) — R{p) relationship could still be 
observed on average. For these cases alone, the estimated 
temperatures do not deviate more than 10% as shown in 
Table IIIIl However, when the resolution is worse than 
5 sites uncontrolled systematic errors dominate and the 
scheme fails. 
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FIG. 8. The fluctuation-dissipation thermometry scheme remains apphcable to in-situ density experiments which have a 
resolution of a few sites. The system consists of 3D optical lattice with A^ = 125, 000 bosonic Rb atoms. Increasing the 
binwidth from 1.0 to 5.0 in steps of 1.0 increases the systematic error in the dissipation term L{p), but the temperature 
estimate remains realiable. Top panels show the dissipation term L{p) and bottom panels show the L{p) vs R{p) curves. The 
parameters in the left column are U/t = 10, T/t = 1; and in the right column U/t — 50, T/t — 1, with a window size ^ = 3. 
We average over a 1000 independent measurements obtained from a QMC simulation. 





FIG. 9. The fluctuation dissipation thermometry scheme can 
still give a reliable estimate for the temperature when it has 
a spread AT/T ~ 20%. The parameters for the 3D bosonic 
*''Rb optical lattice system are N = 125, 000, U/t = 10, T/t = 
1, ^ = 3, and we took 20 independent measurements from a 
QMC simulation. 



FIG. 10. The fluctuation-dissipation thermometry scheme 
remains valid for a 3D bosonic *^Rb system with total particle 
number A^ = 125, 000 and spread AN of about 2,000. The 
parameters are U/t = 10, T/t = 1, ^ = 3, and we took 20 
independent measurements from a QMC simulation. 



n 


-site resolution 


estimated 




(binwidth) 


temperature [t] 




{U/t=10 


T/t=iy. 




1 


0.977 ± 0.007 




2 


0.990 ± 0.006 




3 


0.997 ± 0.006 




4 


1.016 ± 0.007 




5 


1.048 ± 0.007 




( U/t=50 


T/t=iy. 




1 


0.994 ± 0.008 




2 


1.014 ± 0.007 




3 


1.032 ± 0.007 




4 


1.045 ± 0.008 




5 


1.094 ± 0.008 



Appendix C: Fluctuations in temperature T and 
total particle number A'^ 



TABLE III. The estimated temperatures obtained in simu- 
lations for various resolutions. We take a 3D optical lattice 
system with A'^ = 125, 000 ^'^Rb atoms and a window size 
^ = 3. We average over 1000 independent measurements ob- 
tained from a QMC simulation. In the top panel, U/t = 10 
, T/t = 1; in the bottom panel, U/t = 50, T/t = 1. The es- 
timated temperatures do not deviate more than 10% despite 
the systematic error involved. 



Different runs of an experiment will have small tem- 
perature and particle number changes. We model this by 
averaging over simuations where the temperature T and 
the total particle number N fluctuate. First, by letting T 
fluctuate within a spread AT/T of about 20%, wc observe 
no qualititive difference in the thermometry scheme, as il- 
lustrated in Fig.|n]for the case of U/t = 10, T/t = 1 where 
we find an estimated temperature T/t — 0.973(47). 



Second, we let the chemical potential fi fluctuate such 
that the total particle number N fluctuates with a spread 
AA^. We observe that a spread up to 1-2% for the 
case of U/t = 10, T/t = 1,N = 125,000 stifl yields 
a quantitatively reliable estimate for the temperature, 
T/t ~ 1.040(26). However, this scheme will have system- 
atic errors for fluctuations bigger than 3%. 
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